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Abstract. The equilibrium equations for the isotropic Kirchhoff rod are known to form 
an integrable system. It is also known that the effects of extensibility and shearability 
of the rod do not break the integrable structure. Nor, as we have shown in a previous 
paper does the effect of a magnetic field on a conducting rod. Here we show, by means 
of Mel'nikov analysis, that, remarkably, the combined effects do destroy integrability; 
that is, the governing equations for an extensible current-carrying rod in a uniform 
magnetic field are nonintegrable. This result has implications for possible configurations 
of electrodynamic space tethers and may be relevant for electromechanical devices. 

AMS classification scheme numbers: 74K10, 78A30, 74H65 
1. Introduction 

The geometrically exact static equilibrium equations for a uniform symmetric (i.e., 
transversely isotropic) elastic rod are well-known to be completely integrable pQ. In fact, 
there is a close relationship between these equations and those describing the dynamics of 
spinning tops. It is also known that some perturbations of the rod equations are integrable, 
but that others are not. For instance, anisotropy of the cross-section [2] and intrinsic 
curvature [3] destroy integrability, as does the effect of gravity [4] , but extensibility and 
shearability [5] do not, nor does the effect of an external force due to a uniform magnetic 
field @. 

In this paper we show by means of a perturbation analysis that the combined effect of 
extensibility and magnetic field, remarkably, leads to nonintegrability, even though each 
of these effects individually would not destroy integrability. The results may be relevant 
for the study of (localised) spatial configurations of electrodynamic space tethers, i.e., 
conducting cables that exploit the earth's magnetic field to generate thrust and drag 
(Lorentz) forces for manoeuvring [3 [8] . 

The perturbation theory we use is the Hamiltonian version of Mel'nikov theory as 
developed by Holmes & Marsden [9]. This theory considers the break-up of a homoclinic 
orbit of the unperturbed integrable system by detecting transverse intersections of stable 
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and unstable manifolds of a perturbed saddle-type solution. Specifically, the so-called 
Mel'nikov integral measures the distance between stable and unstable manifolds. Simple 
zeroes of this integral correspond to transverse intersections which imply complicated 
horseshoe dynamics and hence nonintegrability. 

From standard results of dynamical systems theory one would also generically expect 
in this situation the existence of a multiplicity of multipulse homoclinic orbits [TO] , 
corresponding to arbitrary many localised solutions of the rod, and we show numerically 
that this is indeed the case. 

Mel'nikov analyses have been applied before to prove nonintegrability of anisotropic 
rods [2], intrinsically curved rods [JT] and heavy rods [4|. The first two of these studies 
employed a formulation in terms of Deprit-Andoyer variables. Here we use a formulation 
in terms of the more common Euler angles. Multimodal configurations for anisotropic and 
intrinsically curved rods have been investigated numerically in [T2], [13] and [3], respectively. 

To apply perturbation theory in the case that two effects are simultaneously present one 
has to make assumptions on the relative scale of the two effects. Since Mel'nikov analysis 
requires closed-form expressions for the homoclinic orbit and these expressions are more 
readily obtained for the extensible rod than for the magnetic rod, the Mel'nikov theory 
is applied with the magnetic field as a perturbation of the extensible rod. However, we 
present numerical evidence, in the form of chaotic Poincare sections, that suggests that 
integrability is also broken in the opposite scaling, i.e., if extensibility is viewed as a 
perturbation of the magnetic rod. 

The structure of the paper is as follows. In Section [2] the governing equations are 
presented as a non-canonical Hamiltonian system. In Section [3] Hamiltonian Mel'nikov 
theory is briefly reviewed. In Section @] the system of equations is reduced to a canonical 
system by using the Casimirs of the Poisson bracket. In Section [5] the homoclinic orbits of 
the unperturbed system are calculated, after which the Mel'nikov analysis is performed 
in Section [61 Multimodal homoclinic orbits and fractal Poincare plots, the signatures of 
spatial chaos, are computed in Section [7] and shown to persist in regions of the parameter 
space well away from the asymptotic region where the Mel'nikov result is valid. Section [8] 
closes this study with some concluding remarks. 

2. The Cosserat Theory of Elastic Rods 

2.1. Kinematic equations 

In Cosserat theory a rod is characterised by a space curve r(s), describing the 
centreline of the rod, and an attached right-handed orthonormal triad of directors 
{di(s), d 2 (s), d 3 (s)}, describing the varying orientation of the cross-section [U]. Here 
s is an arbitrary parameter. 

On introducing a right-handed orthonormal frame {ex, e.2,e^\ fixed in space we can 
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write 

di = Re h (2.1) 

where R is a rotation matrix, i.e., an element of the Lie group S0(3). It is convenient to 
introduce the 'hat map' isomorphism of the corresponding Lie algebra so (3): 

' — a 3 a 2 ' 

1R 3 — > so (3) : a = (ai, a 2 , a 3 ) i— > a = a 3 — aj , (2.2) 

l — a 2 ai J 

so that db = a x b. Differentiating ( 12. ID then gives 

<|J = = i2'i2 T di =: iidi = u x d h (2.3) 

where the prime denotes differentiation with respect to s and n is the vector of generalised 
strains associated with bending and twisting. The body components Ui = u ■ d^ are the 
curvatures (i = 1, 2) and the twist (z = 3) of the rod, for which, from (12.31) . we can write 

Ui = ^ijk d 'j ■ d k- (2-4) 

The second vector of strains v is given by 

r' = v. (2.5) 

The body components v% = v ■ d%, V2 = v ■ c?2 and v 3 = v • ^3 are the strains associated 
with stretching and shear. For an unshearable rod we have v\ = = v 2 . For an inextensible 
rod we have \r'\ = 1. The strain t> 3 actually represents the ratio of deformed to reference 
volume. Since a real rod cannot be compressed to a point, it is natural to impose the 
condition v 3 > 0. For an inextensible and unshearable rod the centreline equation (12.51) 
becomes 

r' = d 3 (2.6) 

and the parameter s can be interpreted as the arclength of the rod. 

It will be convenient in the following sections to express components of vectors with 
respect to the director (or body) frame; for any vector p the triple of components 
(p ■ di,p ■ d 2 ,p ■ d 3 ) will be denoted by the sans-serif symbol p. 

2.2. Constitutive relations 

We assume the rod to be hyperelastic, i.e., we assume that there is a strain energy 
density function W = W (u — u , v — v , s) such that the components of the force 
n = (m, 712, n 3 ) and moment m = (mi, m 2 , m$) in the body are given by 

dW J dW 
rrii = — — and n { = — — . (2.7) 

OUi OVi 

Here uq and vq describe the configuration of the unstressed rod. 
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We will consider the important special case, often called linearly elastic, where the 
strain energy is quadratic in the strains: 

W (u, v) = \b x u\ + \b 2 u\ + \cu\ + \ev\ + \ Jv 2 2 + \k («s - l) 2 , (2.8) 

where B\ and B 2 are the principal bending stiffnesses about d\ and d 2 , respectively, and 
C is the torsional stiffness about c£ 3 . The constants H and J are the transverse shear 
stiffnesses and K is the axial stiffness. In the case of an isotropic rod, B\ = B 2 =: B and 
H — J. 



2.3. Equilibrium equations 



The equilibrium equations for the internal force n and moment m in an elastic rod are 



[14j 



n' + f = 0, 

r' x n + l 



m 



0. 



(2.9) 
(2-10) 



where we have allowed for external distributed loads / and I. The only distributed load 
we shall consider is that due to a magnetic field, in which case / is given by the Lorentz 
force, while 1 = 0. Assuming the rod to be conducting and to carry a current I = Ir' the 
Lorentz force experienced when placed in a magnetic field B is 

f = IxB = Ir'xB. (2.11) 

We assume both current and magnetic field to be uniform and align the e 3 vector of the 
fixed frame with the field, so that B = Be^. Let A = IB, then the equilibrium equations 
when written in the director frame take the form of a non-canonical Hamiltonian system: 



(2.12) 



/ 


m 








f 


m 


n 


e 3 


\ 




n 




= J (m, n, e 3 ; A) VW(m,n) 


with J = -J T = 




n 


Ae 3 







V 


e 3 


J 






\ 


e 3 








/ 



with Hamiltonian 

H (m, n) 



1 

-m 

2 



u(m) + -n 



v n 



;2.13) 



where d 3 = (0, 0, 1). In writing down this equation we have used (12.31) . (12. 5ft and the hat 
isomorphism (12. 2p . This formulation follows that of j6] but now allows for the effects of 
extensibility and shearability. As in the inextensible/unshearable case, the Hamiltonian 
is the same as that of the non- magnetic rod (cf. [13]): the effect of the magnetic field is 
only present in the structure matrix J . 
By introducing the Poisson bracket 



{/,<?} 



(m,n,e3) 



-m • (V m / x V m ^) - n • (V m / x V n g + V n / x V m ^) 

- e 3 ■ (V m / x V e:j g + V e3 / x V m g) - Ae 3 ■ (V n / x V n < 



;2.14) 



evolution of field 



effect of field 
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on (m, 11,63) the equilibrium equations can also be written as 

m' = {m, H} (m n e3) = m x u + n x (d 3 + v) , (2.15) 

n ' = ( n > ^}(m,n,e 3 ) = n x u + Ae 3 x (d 3 + v) , (2.16) 

e 3 ={e 3l % inie3) = e 3 xu. (2.17) 

Note that the Poisson bracket is an extension of the usual Kirchhoff bracket by two 
terms. The first is a semidirect term [16] describing the evolution of the magnetic field 
in the director frame (this term does not affect the force and moment balance since 
the Hamiltonian is independent of e 3 ). The second is a cocycle known as a Leibniz 
extension [TTj describing the effect of the magnetic field on the conducting rod. Note 
that (I2.17P is just the equation e' 3 = written in the director frame. 

Note on notation: an ornamented gradient symbol is used for gradients with respect to 
the indicated fields, while the unornamented symbol will always denote the gradient with 
respect to all three fields (m, n, e 3 ). 

The non-canonical system (12.15p - (l2.17j) has three Casimirs, given by 

Ci = -n-n + Anve3, (2.18) 

C 2 = e 3 • n, (2.19) 

C 3 = e 3 ■ e 3 . (2.20) 

The magnitude of force is not conserved if a magnetic field is present, but as a result 
of rotational symmetry the force component in the direction of the field is conserved. 
Casimir d 2 . 1 8 [) does not seem to have a physical interpretation. Naturally, C 3 = 1. 

In the special case of inextensible/unshearable and isotropic rods two first integrals 
emerge: 

h = Bm-d 3 (if J = H and B X = B 2 = B), (2.21) 
J 2 = n • m + B\e 3 • d 3 (if J = H = K = and B X = B 2 = 5)(2.22) 

As in the Kirchhoff case, the first of these expresses conservation of twist in the rod. 
The second integral does not seem to have a physical interpretation, but in the limit 
A — > it reduces to the familiar conservation of torque about the loading axis. It is 
a straightforward task to check that all the first integrals (12.13j) . (12.211) and ( 12.221) are 
independent and in involution with respect to the Poisson bracket (I2.14p . and therefore 
the system is completely integrable. 

Note that the twist integral I\ only requires isotropy, while the integral I 2 requires 
isotropy and inextensibility/unshearability. Indeed, it is the main purpose of this paper 
to show that the combined effect of inextensibility/unshearability and magnetic field leads 
to nonintegrability. For this we use Mel'nikov theory, which we review next. 
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3. Mel'nikov theory 

For Hamiltonian systems Mel'nikov's original perturbation analysis needs to be 
adapted. We use the results for two-degrees-of- freedom systems presented in [HIE]. 
Consider a Hamiltonian that depends on a small parameter e in the form 

H (q, p, if, I) = Ho (q, p, I) + eHx (q, p, <p,I) + (e 2 ) , (3.1) 

where (q,p) are conjugate variables and (J, <p) are action-angle variables such that Tii is 
27r-periodic in ip. For e = ip is a cyclic variable (hence / is a first integral) and Hamilton's 
equations are completely integrable. We assume that the unperturbed Hamiltonian TCo 
satisfies the following two conditions: 

(i) For some / = I Hamilton's equations corresponding to 7{q possess a homoclinic orbit 
(<? (t) iV CO) to a hyperbolic fixed point (qo,po) at Hamiltonian level h = Ho(q,p, Io)- 

(ii) The frequency 

:= j^-(q,PJo) (3.2) 

of the unperturbed system satisfies |co>o| > v > for some v e K. and Vt G (— oo, +oo). 
This condition means that ip is a time-like variable and allows the unperturbed system 
to be reduced to the (q,p) space with ip as the independent variable. 

Now define the Mel'nikov function 

MM = f °° {h ,—\ dt, (3.3) 



LO 



where the canonical Poisson bracket {/, g}{ q , p ) = q~q~ ~ f^ff ^ s evaluated at the 
homoclinic orbit and ip(t) = f*u> (t) di + <p . We then have the following result: 

Theorem 3.1. For e ^ sufficiently small, if M. (<^o) has simple zeroes <po, that is, 

dM 

M(f ) = and — (<p )^Q, (3.4) 
o<p 

then the stable and unstable manifolds of the perturbed hyperbolic invariant set (a periodic 
solution in the four-dimensional system) intersect transver sally. If, on the other hand, 
A4 (ipo) is bounded away from zero, then the manifolds do not intersect. 

Proof o/IO See [21 El [18] . □ 

The existence of transverse homoclinic orbits implies that the 'dynamics' near the 
hyperbolic saddle is 'chaotic' in the sense that the following holds: 

Corollary 3.2. The Poincare map associated with Tt on the homoclinic level set 7i _1 (/i) 
has a hyperbolic, non-wandering Cantor set on which the map is conjugate to a Bernoulli 
shift of finite type. 

Proof o/E3 See [2] and [18]. □ 
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This in turn implies 

Corollary 3.3. The Hamiltonian 7i has no analytic conserved quantities independent of 
7i itself, i.e., the corresponding Hamiltonian system is nonintegrable. 



Proof of [Q See [2j. 

4. Reduction of the magnetic rod equations to a canonical system 



□ 



In this section the three Casimirs (I2.18l) - (l2.20p are used to reduce the nine- 
dimensional non-canonical Hamiltonian system (I2.15I) - (I2.17I) in (m,n,e3) to a six- 
dimensional canonical Hamiltonian system in terms of Euler angles and their canonical 
momenta (q,p) = (#,?/;, 0,p6>,Jty'£V>)- The reduction follows [6j but now allows for 
extensibility and shearability of the rod. The reduction in [5] was shown to be canonical on 
the condition that the force n and magnetic field B are not aligned. This result trivially 
extends to the present case as the effects of inextensibility and shearability do not enter 
the structure matrix, only the Hamiltonian. 

Let 

^ cos 9 cos cos ip — sin sin ip cos 9 cos sin ip + cos ip sin — sin 9 cos 6 \ 
— cos 9 sin cos ip — cos sin ip — cos 9 sin sin ip + cos cos ip sin 9 sin 
sin 9 cos ip sin 9 sin ip cos 9 

be a parametrisation of the rotation matrix R in (12 .ip in terms of Euler angles. Here 9 is 
the angle the tangent to the rod makes with the magnetic field, ip is the azimuthal angle 
about the field direction and is the twist angle about the centreline of the rod. It follows 
that for the triple e3 we have 



R 



\ 



J 



e3 (g) = R (q) k = (— sin 9 cos 0, sin 9 sin 0, cos ( 



(4.1) 



where k = (0, 0, 1) T . On inserting the Euler angles into the strains (12.41) and (I2.5P and 
using the constitutive relations (12.71) the moments are found to be 



m 



/ 


mi 


\ 




f 




9' sin — ip' sin 9 cos 


)) 






m 2 








B 2 ( 


9' cos + ip' sin 9 sin 


) 


= Lp, 


V 


m 3 


J 




V 




C(0' + ?//cosfl) 


) 





(4.2) 



where 



sin 9 



' sin 9 sin < 
sin 9 cos < 

V 



- COS i 

sin 




cos 9 cos 
— cos 9 sin i 
sin 9 



\ 
J 



and the canonical momenta defined by pg = dW(q,q')/d9', p^ = dW(q,q')/dip' , p^ = 
dW(q,q')/d(p', with W(q, q') = W(u(g, q')) in terms of the strain energy function W 
defined in (12131) . 
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For the force we can write n = R (q) w (q,p), for some non-constant triple w. By 
decomposing w into parts perpendicular and parallel to k and using the Casimirs (12.181) 
and (12TT9]) we obtain [6] 



C, 



V 



- sin 9 cos i 
sin 9 sin 
cos 9 



\ 



/ 



+ J2C l - Ci - 2\p, 



^ cos 9 cos cos if) — sin sin if) ^ 
— cos 9 sin cos if) — cos sin ip 
sin 9 cos if) 



(.4.3) 



It was shown in [6] that the transformation (I4.ip - (I4.3I) is canonical provided that the 
force and magnetic field (or 63) are not aligned, i.e., 

2C l - C\ - 2X Pi> ^ 0. (4.4) 

(It was also shown that if (14.41) holds anywhere along the rod it holds everywhere, and in 
that case all solutions are straight twisted rods, aligned with the magnetic field.) In the 
linearly elastic case the Hamiltonian (I2.13P transforms into 



H(9,if),(f),p e ,p i} ,p ( j ) ) 



2B 1 



+ 



ml 



rr. 



2B 2 + 2C + 2H + 2 J + 2K + ^ 



(4.5) 



with rrii, given in terms of the canonical variables by (I4.2p and (|4.3|) . In the isotropic 
case (Bi = B 2 =: B , H = J) the Hamiltonian reduces further to 

2 



li{0,il),pe,p^,p^) 



1 



+ 



2B 
+ ( Co 



2B 
1 

~K 



Pi> - p<f> cos 9 



sin 9 



C2 
2 



1 

K 



cos 9 



+ C 2 cos 9 

cos 9 + 1 J sin 9 cos if> \J 2C 1 — C 2 — 2\p^ 

+ \{~K~l) sin2 9 cos2 ^ {2Cl ~ °* " 2Ap ^ ~ 7 Pv " (46) 



where we have dropped the p 2 ^ term, which is constant since is a cyclic variable: 
p^ = m 3 = Ii/B. If, in addition, the rod is inextensible and unshearable (H = J = K = 0) 
then Hamilton's equations corresponding to (14. 6ft have I 2 in (I2.22f) as a first integral, which 
in canonical variables takes the form 

'p (j> -p^cos9\\ 



I 2 = \B cos 9 + C 2 p^ — \ 2C l — C| — 2\p^ I p g sin if) — cos if> 



rendering the system completely integrable. 

Finally, we use the constants C 2 and m 3 
quantities by setting 



h/B 



sin 9 

Ptj, to introduce dimensionless 



m 3 _ p e 
t = s—, p e = — , 
B m 3 



Pi> 



m 3 



7 = C 2 



1 

K 



J 



m 



Am 3 
= C| ' 

m 3 
\J BC 2 



2d 



c 2 2 



c| 



(4.8) 
(4.9) 



so that the dimensionless Hamiltonian H = 7iB/m\ and integral I 2 = I 2 /{C 2 m 3 ) become 
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- 1_ 2 1 fp^-cos9\ cos6 7COS 2 ^ 



2™ 2 V sin / m 2 2m 2 
+ -^7(7 cos6» + 1) smOcostpJu- 2\pJ (4.10) 

+ 77^7 sin 2 6 1 cos 2 ip(fJ, — 2Xp^) - -^p<d> 
A cos 9 I - _ 1 ( _ . , , / 1 — pit co$9\ \ . . 

^ + - ~ 2A ^ l^ sin ^ ~ cos ^ I — — ) ) ■ t 4 - 11 ) 

Remark 4.1. J2 zs not on/?/ a first integral of the canonical equations generated by Ti for 
7 = 5 = 0. It is also a first integral for the case A = 0. This is seen from the more familiar 
form in ( QTU| ) more readily than from the form in (4-H ) obtained after going through the 



and 



reduction. In the absence of a magnetic field the parameter \x is artificial, the result of our 
choice o/e 3 . 

Since it seems impossible to use the integral I2 to reduce the canonical system further 
to a single-degree-of-freedom one, homoclinic orbits are not easily obtained in the general 
case. However, in the non-magnetic case (A = 0) a further reduction is possible and 
homoclinic orbits can be obtained explicitly. This is the topic of the next section. 

5. Homoclinic solutions of the extensible rod in zero magnetic field 

In the absence of a magnetic field (A = 0) the force n is a constant vector, by (12.91) . 
Provided this vector is not zero, we can choose the fixed frame vector e 3 in the direction 
of n (this gives the usual physical meaning to the Euler angles 9, ip, 4>)- According to 
(I2.18P and (12. 19)) the Casimirs satisfy 2C\ = C\ 7^ 0. Hence \x = and the Hamiltonian 
(14.61) becomes 

- _ 1_ 2 1 fp i ,-cos9\ 2 cos9 7COS 2 # 

We are interested in homoclinic orbits, so we assume the rod to be loaded by an end 
force and end moment applied axially to the rod, which is aligned with e 3 as t — > ±00. 
Thus n 3 and m 3 are the end loads and we have p^ = n • xx\jG% = J2/C2 = ^3, hence 
p,p — 1. Hamilton's equations corresponding to (15.1 ft in this case read 

^1 — cos#) 2 (7 cos 9 + 1) sin# 
sin 3 9 m? 

where we have used an overdot to denote differentiation with respect to t. This system of 
equations agrees with that derived in [5] (see also the planar reduction in [T9]). 

The trivial fixed point 9 = of (15.21) corresponds to a straight twisted rod. Non-trivial 
fixed points solve the cubic 

(7COS0+ 1)(1 + cos#) 2 = m 2 (5.3) 



9 = p e and p e = - ^ .77 — V , (5-2) 

1 f-l rn ^ 
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and correspond to helical solutions. They exist for 

0<m 2 <4(l+ 7 ), (5.4) 



where the upper limit corresponds to the critical load m c = 2-y/l + 7' for torsional buckling 
described by a pitchfork bifurcation [3]. 

For parameters satisfying (15 .4p the trivial fixed point is a hyperbolic saddle from which 
a symmetric pair of homoclinic orbits emanates. To find these we integrate (I5.2p once to 
obtain 

lp + V{0) = h, ^) = - 1 " C °1 + ^ + ^cos^ (5.5) 
2 K J y J 2 1 + cos# m 2 2m 2 v ; 

with h the 'energy' level. The corresponding phase portrait is shown in Figure 1(a) (while 

Figure 1(b) shows a typical phase portrait for m > m c ). Setting u = cos 6 and noting that 

the energy of the homoclinic orbits is 

m z V 2 

we can solve (15.51) for u and integrate to get 
m r® du 

t = — I ; (5.6) 

y/l Ju(o) (l-u)y/g (u) 

where 

g (u)=u 2 + 2u(l + -) +1 + --—. (5.7) 

V 7/ 7 7 

The quadratic g (u) has roots 



u± = -(l + -] ±-V^+ 



7m 2 



Substituting the upper and lower bounds for m in (15. 4p into these roots gives 

3 + — J < w_ < - ( 1 + - J and - 1 < u+ < 1. (5. 



2 7/ V 1, 

Thus, the roots of g are always distinct, i.e., w_ 7^ u + . The integral (15.61) and the limits 
of integration can be simplified as 

m f u{t) du . . 

t=—^ I . =;. (5.9) 

V7 Ju+ {l — v)\J(u — «_) (u - u + ) 
The substitution u = U- + (u+ — uJ) cosh 2 z turns the integral into 

t = r / 5- with fc = -± . (5.10) 

(1 — m_) v /7 Jo 1 — A; cosh 2 1 — tt- 

Note that since m_ < — 1 and — 1 < u + < 1 we have that > 1 for all parameter values. 
Thus the integral (I5.10p is never singular and may be solved by using the identity 

dz 



1 — k cosh z 




(5.11) 
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Figure 1. (a,b) Phase portraits of (|5.5p for 7 = 1, m = 1.7 < m c and m — 2.2 > m c , 
respectively. (c,d,e) 0, and ip as a function of scaled arclength f for the homoclinic 
orbit of the extensible (solid) and inextensible (dashed) rods at 7 = 1, m = 1.7. 



Hence, the homoclinic solutions are given by 



cos6 = u^ + («+ - M _)cosh 2 ^tanh" 1 ^^^^tan ^ ^zj^g lJl j 

p e = 0. (5.13) 

This solution agrees with that derived in [51 Appendix], where it is expressed in terms of 
a natural logarithm rather than hyperbolic functions. In the limit of small extensibility, 
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i.e., 7 — > 0, (I5.12p recovers the expression [201 Eq. 69]: 

9 = cos' 1 (u + (1 - M ) tanh 2 ^ ^ tj J and p e = 9, (5.14) 



where u + — > u = m 2 /2 — 1, u_ — > — oo and fc — > 1 as 7 — > 0. Figures l(c)|l(d)|l(c' 
compare both homoclinic orbits. 

The derivative of the angle ■?/> is given by 

i> = J-— n- (5-15) 

1 + COS 9 



Figures l(c)|l(d)|l(e) show plots of 9, pe and ip, which we will need in the Mel'nikov 



analysis. 

6. Mel'nikov theory applied to the magnetically perturbed extensible rod 

In order to express the Hamiltonian (14.101) in the form (13. II) for use in the Mel'nikov 
analysis we introduce a small parameter e and write 

/i = ae 2 and A = be 2 , (6.1) 

where a and b are positive and O (1), as are the other parameters m, 7 and 5. Then the 
Hamiltonian takes the form 

H (9,i/j,pe,p i ,) = H (9,p e ,p^) + eU x (9,ij;,p9,p^) + O (e 2 ) , 

where the unperturbed Hamiltonian 7i is given by (15. ip and the first-order perturbation 
is given by 

(9,ip,p e ,Pih) = -At (7COS6 1 + 1) sin6 l cos^v /a _ 2bpJ. 
For the frequency at the homoclinic orbit (15.121) we find 

dn 



°Pl> hom 1 + COS^ 

Since 9 in the homoclinic orbit is bounded away from tt, this ujq is well-defined and bounded 
away from zero. The two conditions (i) and (ii) in Section [3] are therefore satisfied and 
Mel'nikov theory can be applied. 
The required partial derivatives are 

dH . n ( 1 (l+7cos( 



09 V(l + cos^) 2 m 2 

dHo _ 
ope 

8ujq sin 9 

~d9 ~ (1 + cos#) 2 ' 



dp e 



0. 
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Figure 2. Plots of the (normalised) Mel'nikov integral at m = 1.7 and 7 = 1, 1.2, 1.4, 
1.6, 1.8 and 2 (from top to bottom at ipo — it/2). The integral is identically zero for a 
value of 7 somewhere between 1.8 and 2. 



dl~~L 1 

= — o ( cos # + 7 cos 26) cosip\/a - 2bp ll ,\ 
06 m 2 

^ = 0. 
On using (13.31) and the identity 

^o,^-} = — {^o,^i} (e ,p fl ) - -4 {^o,^o} (e .p e) , (6.3) 

and writing if) (t) = ip (t) + V>o> with such that ip(0) = 0, the Mel'nikov integral is found 
to be 

— y/a - 2b' f +oc _ 

M(ip Q ) = 5 sin-^o / pgsmip [(1 + cos6>) (cos# + 7 cos 20) (6.4) 

i-00 

+ sin 2 0(1 + 7 cos 6)] dt. 

Here we have dropped the cos ipo term which by symmetry does not contribute since 6 is 
an even function of t while pe and ip are odd functions of t (cf. Figured]). 

Generically the Mel'nikov integral will have simple zeroes provided a — 2bp^ 7^ 0. This 
condition is no restriction as it corresponds to the non-alignment condition (14. Ah which 
is assumed throughout. There may be special parameter values for which the Mel'nikov 
integral is zero. In these exceptional cases the system will remain nonintegrable but the 
intersection of the stable and unstable manifolds will be nontransverse. 

A plot of the Mel'nikov integral (16.41) is shown in Figure [2] for various values of 7 
confirming the existence of simple zeroes. We note that for 7 close to 1.8 the integral 
is identically zero, indicating a nontransverse intersection of the stable and unstable 
manifolds. 
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7. Numerical results 

By Corollary (13.21) we expect the equations for the extensible magnetic rod to be 
chaotic, i.e., to contain a horseshoe (at least for small A). To confirm this we present 
Poincare sections in Figure [31 Each panel in the figure shows one orbit with starting values 
6 = 0.1, pe = 0.5, ip = at fixed Hamiltonian level h = 0.9. Solutions were computed 
using the 8th-order Dormand-Prince code D0P853 [2T] with relative error tolerance set 
to 10~ 12 and intersections were recorded with plane of section given by sinip = 0. Plots 
represent 10000 Poincare iterates taking runs up to t = 84000 in which the Hamiltonian 
was found to be preserved to within 3.6 x 10~ 9 . Clearly visible is the break-up of regular 
closed orbits into the typical fractal sets of chaotic systems as A is varied. 

Given the existence of a transverse homoclinic orbit in a chaotic system one expects 
from standard results from dynamical system theory also the presence of higher-order 
(multipulse) homoclinic orbits that correspond to solutions that pass near the saddle 
solution (and the unperturbed homoclinic orbit) multiple times before closing up at 
the saddle \±0\. This was found to be the case for anisotropic rods in [13] . Multipulse 
homoclinic orbits for the magnetic rod are displayed in figure HI They were obtained by 
means of the shooting method discussed in [13] . 

8. Discussion 

We have shown that the equilibrium equations for an extensible and shearable 
conducting rod in a uniform magnetic field are nonintegrable and exhibit chaotic solutions. 
This is surprising as the effects of extensibility/shearable and magnetic field individually 
do not lead to a breakdown of integrability of the classical isotropic rod equations. Our 
system is unlikely to be unique in showing this lack of 'additivity' of integrability, but we 
have not seen it reported of other systems before, physical or otherwise. 

To prove nonintegrability of the system we used Mel'nikov's method with the equations 
scaled in such a way that the unperturbed rod is extensible/shearable but non-magnetic 
and the magnetic effect forms the perturbation. This is necessary because the Hamiltonian 
system for a (inextensible) magnetic rod, although integrable, cannot be explicitly reduced 
to a one-degree-of-freedom system (at least not globally in terms of Euler angles) and 
therefore the required explicit expressions of the homoclinic orbit are not available. It 
is reasonable to expect, however, that integrability is broken more widely in parameter 
space, including regions of small extensibility (7 <C 1) and large magnetic field (A = 0(1)). 
Figure [5] gives numerical evidence for this in the form of chaotic Poincare plots computed 
for different parameter values. Each panel in the figure shows one orbit with starting 
values 9 = 0.1, pg = 1, ip = at fixed Hamiltonian level h = 0.8. Plots represent 10000 
Poincare iterates taking runs up to t = 81000 in which the Hamiltonian was found to be 
preserved to within 3.9 x 10 -9 . Nonintegrability for these parameters is also confirmed by 
the multipulse homoclinic orbits for small 7 and 8 shown in Figure El 
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Figure 3. Poincare sections for sinip = at energy level h = 0.9 for varying (small) 
values of A. (m = 1.7, fi — 0.4, 7 = 3, 6 = 3.) 




-1 -0.5 0.5 1 



Figure 4. Two-pulse (solid) and four-pulse (dashed) homoclinic solutions for small A 
and fj, (A = 0.001, n = 0.002, m = 1.7, 7 = 1, <5 = 1). 



As the chaotic solutions correspond to spatially complex configurations of the rod, our 
results may be relevant for electrodynamic space tethers [7J |8] and, at an entirely different 
scale, for beams or ribbons as part of micro- or nanoelectromechanical devices such as 
sensors, resonators, inductors and actuators [22]. For instance, there is significant interest 
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Figure 5. Poincare sections for sin^ = at energy level h = 0.8 for varying (small) 
values of 7. (to = 1.7, S = 0.25, A = 2, fj, = 4.) 



1.6 
1.4 
1.2 
1 

<as 0.8 
0.6 
0.4 
0.2 





-0.5 




t 



0.5 



Figure 6. Two-pulse (solid) and three-pulse (dashed) homoclinic solutions for small 7 
and S (7 = 0.01, S = 0.01, m = 1.7, A = 0.5, (j, = 1). 



in nanosprings of small pitch because they allow for large magnetic flux densities [23] . 

For sufficiently slender elastic structures localised (i.e., homoclinic) solutions are the 
preferred mode of deformation [2D]- We have presented preliminary numerical results 
showing that in addition to the transverse homoclinic orbit guaranteed to exist by 
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Mel'nikov theory there exist multipulse homoclinic solutions. It would be interesting to 
study the bifurcation behaviour of these localised solutions as physical parameters are 
varied. We intend to take this up in a future publication. 
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